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We present a new method to measure the guided, radiated and total decay rates in uniform waveguides. It is 
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It is well-known that the Spontaneous Decay Rates 
(SDRs) of emitting sources such as atoms or quantum 
wells can be largely modified in optical micro-cavities. 1 
Controlled spontaneous-emission plays a key role in a 
new generation of micro and nano-optical devices. High- 
performance micro-cavity lasers, for instance, have been 
already experimentally demonstrated. 2,3 SDRs modifi- 
cations are also expected to impact the performance of 
optical waveguide amplifiers as the guided mode area 
radii of these devices become smaller 4,5 than s» l/im. In 
order to characterize one dimensional nano- and micro- 
optical devices one key decomposition of the total decay 
rate (tq 1 ) into two components was introduced in Rcf. 
[6]. The decay rate was divided into the guided modes 
(t^ 1 ) and into the radiated modes (t" 1 ), where the to- 



tal decay rate is r 1 = t q 



In contrast to large 



devices, where just the total decay rate must be consid- 
ered for their characterization, the modeling of nano- and 
micro-optical devices requires the measure of both com- 
ponents of the decay rate. The usual way to determine 
the SDR of an emitting source embedded in a uniform 
waveguide consists in the measure of the exponential de- 
cay rate of the Amplified Spontaneous Emission (ASE) 
output power when the pump source is switched off. 7 
The SDR is given by the exponential decay coefficient of 
the ASE output power. 

Two natural questions arise here:(l) How can r g and 
r r be measured? (2) What is actually measured when 
using the classical method to determine the SDR? In 
this Letter we answer these two questions. We show that 
the classical method to measure To in uniform waveg- 
uides gives r r if used in long length waveguides and ac- 
tually To in short waveguides (assuming no reflections at 
the waveguide ends). We also show how these measures 
are modified in lossy mediums, i. e., when a background 
loss coefficient is incorporated into the rate and prop- 
agation equations. Then we show how these ideas are 
useful on devices of practical interests. Three cases are 



considered: Erbium Doped Waveguide and Fiber Ampli- 
fiers (EDWAs and EDFAs) and Semiconductor Optical 
Amplifiers (SOAs). 

Theory: decay rate measures and background loss in- 
fluence. We employ the analytical solution for the longi- 
tudinal z dependence of the rate equations presented in 
Ref. [5] to investigate the measure of the SDR in uniform 
waveguides. The analytical solution presented there is 
valid only for waveguides in which the excited state pop- 
ulation of the emitting source, N2(z), is constant along 
the fiber. Since the measure of the decay rate is per- 
formed when N 2 — > along z, this approximation is 
valid when measuring tq. The rate equation is: 5 



dN 2 (z,t) 
dt 
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where A/"i + A/2 = 1 are the normalized population of the 
upper and lower levels of the emitting source, t is the 
spontaneous lifetime of the upper level, p is the number 
density of active ions, S is the doped region area, and 
a„ and y„ are the absorption and gain constants. The 
propagation equation is: 



dP n (z,t) 
dz 



&n &loss\ 

x P n (z,t) + 2j n AvAf 2 (z,t)}, (2) 



where P n (z,t) is the optical power (in photons per unit 
time) at location z of the nth beam with wavelenght 
centered at A„ (n < M), u n = 1 for forward travel- 
ling beams and u n = — 1 for backward travelling beams, 
&ioss is the attenuation coefficient given by the back- 
ground loss of the fiber glass host, Av is the frequency 
interval between two successive wavelengths considered 
in the model, and the factor 2 in the last term stands for 
two possible polarizations. Solving Eqs. Q and (J2J) for 
N 2 (z,t) = A/2(t), i- e., A/2 constant along z, the output 



1 



power is 

Pn Ut (t) 



PT{t)G n (t) + ZA/^Ai/ [G n (t) - 1] , (3) 



where 



G n (t) = 



7n)-Mi(*) - a„ - a ioss 



The rate equation is 
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where 



P»(i) 



ln[G„(t)] 
HG n (t)} 



[G n (t)-l]+2Af* p Ai 



(7) 



Here P° U *W = Pn(-M) and P™(t) = P n (0,t) are the 
output and input power of the nth beam, G n (t) is the 
linear gain, Af^ p is the spontaneous emission factor for 
the nth mode, and L is the doped fiber length. 

Measuring the decay rate. In the classical method 7 
to determine the SDR, the input power is turned off 
(P" l (t) = 0) and the useful data is collected when the 
concentration of excited ions is low (A/2 <C 1). With these 
two conditions Eqs. © and J2J) become: 



H n (t) 



2 ln AvN 2 {t) 
a n + aioss 

K ut {t) 
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Using Eq. @ we can write Eq. © as 



, (8) 
(9) 
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dt 
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where we introduce the effective background loss coeffi- 
cient for the nth mode (3 n = ai OS s/ (a n + ai OS s). Looking 
at Eq. JHJ) , which is linear in Af 2 (t) , we see that the right 
hand side of Eq. (|10|l is also linear in Af 2 (t). Therefore, 
it can be rewritten as dAf 2 (t)/dt = —Af 2 (t)/T m , where 
r m is what is actually measured by the classical method 
and not To- Since r m is quite cumbersome, we do not ex- 
plictly write it here. But two limiting cases deserve a de- 
tailed study. Case 1: short length waveguides. In this case 
(a n + aioss)L <C 1. With this approximation Eq. IjlUfl re- 
duces to 

dVa(t) _ M 2 {t) 



dt 



(11) 



This result shows that only for short length waveguides 
the classical method 7 furnishes the total SDR of the ion. 
It is interesting to note that Eq. Ijllfl is valid whether or 
not we have background loss (ai OS s 7^ 0). Case 2: long 
length waveguides. Here (a n + ai OS s)L ^> 1. Now Eq. (|10|) 
becomes 



dM 2 (t) 
dt 



M 2 {t) ^M(t) 



(12) 



n=l 



where we have used the decomposition 6 of r 1 in guided 



and radiated modes (r 



r r 1 ) and the fact 



- Y%Li T g}> where r g } = 2^ n Av/pS is the 



that Tg 1 

guided decay rate into the nth mode. For sufficiently low 
background loss (3 n w 0, which implies that the classical 
method now furnishes t^ 1 . 

The fact that r^ 1 in a given mode can be written as 
= 2 r j n Au/pS is pointed out here for the first time. 
It arises naturally from the interpretation given in Ref. 
[5] for this term as being the photons captured by the 
guided modes per unit time. In contrast to Ref. [6], in 
which t~ 1 is given as a function of the dipole moment 
matrix element between the emitting source excited and 
ground states, we use here the easily measurable gain 
constant j n . 

Therefore, the total (tq 1 ) and the radiated (r" 1 ) de- 
cay rates in a given uniform waveguide can, in principle, 
be measured separately employing two different waveg- 
uide lengths, t^ 1 can also be determined by simply sub- 



tracting the later from the former (r 



Anyway, t" 1 is also easily obtained from the waveguide 
intrinsic parameters, as discussed in the previous para- 
graph. 

We observe that in Ref. [5], Section 2.B.3, the results 
above were outlined for ai oss = 0. But now we explicitly 
perform the calculations and reinterpret these results in 
light of the decomposition of the decay rate into guided 
and radiated modes. 

Simulations. We end this letter studying waveguide 
lengths range at which t^ 1 and t" 1 can be measured. 
We use three sets of parameters of practical interest. 
They represent typical EDFAs, EDWAs, and SOAs. It 
is worth mentioning that the term proportional to A/f , 
which is usually included in the modelling of EDWAs 
and SOAs, can be neglected since M 2 <C I. 

EDFAs and EDWAs. In usual EDFAs and EDWAs, 
with optical mode areas larger than « lpm 2 , r^ 1 is 
negligible and the total lifetime is equal to the radiated 
lifetime (tq w t t w 107ns). However, when the optical 
mode area becomes smaller, r" 1 starts to have a non- 
negligible and measurable value. For instance, using the 
parameters shown in Tabled a lnni optical mode ra- 
dius, and Eq. (|1(J|) . we determine the measured lifetime 
as a function of the waveguide length for typical EDFAs 
and EDWAs (assuming perfect detection). Without loss 
of generality, we use an effective gain and absorption con- 
stant along the total transition bandwidth of « 15Thz 
centered at 1.55uto. We observe in Fig. Q that a lpm 
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Tabic 1. Parameters used in the simulations. 
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1.4 x 10 26 



E 



(a) 




40 6U 

length (cm) 

Fig. 1. Simulations of the measured lifetimes as a function of 
the waveguide length for typical EDFAs (a) and EDWAs (b). The 
black lines are for non-lossy waveguides and the light grey ones for 
"loss = 0.3dB/m (EDFAs) 9 and ldB/cm (EDWAs), 8 typical of 
fluorozirconatc EDFAs and silica based EDWAs, respectively. 

optical mode radius is small enough to cause a varia- 
tion of w 5% between the total and the radiated life- 
times, measured at small and large waveguide lengths, 
respectively. Although small, an optical mode radius of 
« 1/iTO is already commercially available in Photonic 
Crystal Fibers. Moreover, it has been shown that waveg- 
uides with high index-contrast (and, consequently, small 
mode areas) have several advantages, 10 which envisages 
the future construction of very-small mode area devices. 
To study how the total and radiated decay rates can be 
measured in such devices, we perform simulations assum- 
ing the values shown in Table \I\for an EDWA, but an 
optical mode area of 0.02pm 2 , which was obtained, for 
instance, in Rcf. [3] . The results are shown in Fig. [21 Of 
course, such highly-confined EDWA would use materials 
which will not necessarily have the parameters shown in 
Table ^ But the graphic in Fig. [2] is illustrative of the 
effects that would always occur at smaller mode areas: 
the difference between the total and the radiated lifetime 
increases (as a consequence of larger t" 1 ), the difference 
between the decay rates in lossy and non-lossy waveg- 
uides increases, and, at last, the waveguide lengths to 
obtain the radiated lifetime with a given accuracy also 
increases. 



SOAs. In SOAs, the condition A/2 + M\ = 1 does not 
hold, but it is easy to show that the method here pre- 
sented to measure the radiated and total lifetime also 
works. We observe that 7„ can be written as cr^Tp /S, 
where is the emission-crosssection at wavelength A n 
and r is the overlap factor between the optical mode 
and the doped region area. 11 Then, t~ 1 is given by 
t- 1 = 2a*T/S. Using a* T = 2.0 x 10" 16 cm 2 for a 
typical SOA, 12 we found that r" 1 becomes compara- 
ble to TjT 1 of « (200ps) _1 only at optical mode radius 
smaller than rs lnm, far from the possibilities of present 
technologies. 

G. Rigolin thanks FAPESP for funding this research. 
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Fig. 2. Theoretical lifetimes as a function of the waveguide lenght 
for an EDWA with 0.02^m 2 optical mode area. Other parameters 
are given in Tab. HI 
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